Method for computation of differential azimuth from spaced-apart gravity component measurements

ABSTRACT

A method is provided for computing the differential azimuth between two sets of spaced-apart gravity component measurement devices in a borehole, employing multiple rotation angles.

BACKGROUND OF THE INVENTION

This invention relates to the computation of differential azimuth between two sets of spaced-apart gravity component measurement devices in a borehole, subject to the constraint that there is no rotation about the borehole axis between the two sets of measuring devices, generally known as accelerometers. Such a structure cannot determine absolute azimuth orientation but only the differential azimuth between the two spaced-apart sets of gravity component measurements. Such apparatus in general is referenced in prior art U.S. Patent Application Ser. No. 06/200,096, Borehole Survey Apparatus and Method, filed Oct. 23, 1980, describing one such a mechanization. That application was processed to the point of allowance of all claims and then abandoned. Another prior art example is described in U.S. Pat. No. 6,480,119, Surveying a Subterranean Borehole Using Accelerometers, issued Nov. 12, 2002. No details are there provided on the development of the equation shown therein for differential azimuth and there are issues with that equation that lead to a conclusion that acceptable results are not consistently obtainable for all possible input conditions.

It is an objective of the present invention to provide a preferred method for obtaining the above cited differential azimuth measurement that is completely described so as to teach and show the method and to provide an improved method verifiable as correct by independent analysis by users. Further, it is an objective to provide two alternative embodiments that are simpler computationally, but demonstrate less geometrical accuracy than the preferred method. These two alternatives are similarly completely described.

SUMMARY OF THE INVENTION

The invention provides a method for computing the differential azimuth between two sets of spaced-apart gravity component measurements or devices, in a borehole, subject to the constraint that there is no net rotation about the borehole axis between the two sets of gravity component measuring devices, the method comprising the steps of:

a) modeling the relationship of the sensitive axes of a first set of said gravity component measurement devices to an earth fixed coordinate set in terms of an inclination rotation angle, for example alpha, from the vertical and an orthogonal highside rotation angle, for example beta, about said borehole longitudinal axis,

b) modeling the relationship of the sensitive axes of a second set of said gravity component measurement devices to said sensitive axes of said first set of gravity measurement devices in terms of a rotation angle, for example tau, a bend direction angle, about said borehole axis, and another rotation angle, for example kappa, a bend magnitude angle, about an axis that is perpendicular to said borehole axis and defined by the bend direction angle tau and a last rotation angle that is the negative of said first rotation angle tau about the said borehole axis,

c) forming equations relating said rotation angles alpha, beta, tau, and kappa to the outputs of said gravity measuring devices,

d) solving said equations to derive the values representative of all of said rotation angles, using measured values of said outputs from said first and second sets of gravity component measuring devices, and

e) computing from said values of said rotation angles the differential azimuth angle between said first set and said second set of gravity measuring devices.

An alternative method includes the steps:

a) modeling the relationship of the sensitive axes of a first set of said gravity component measurement devices to an earth fixed coordinate set in terms of an inclination rotation angle, for example alpha, from the vertical and an orthogonal highside rotation angle, for example beta, about said borehole longitudinal axis,

b) modeling the relationship of the sensitive axes of a second set of said gravity component measurement devices to said sensitive axes of said first set of gravity measurement devices in terms of a rotation angle, for example gamma, about an axis perpendicular to said borehole axis and another rotation angle, for example delta, about another axis that is perpendicular to said borehole axis and also perpendicular to other axis,

c) forming equations relating said rotation angles alpha, beta, gamma and delta to the outputs of said gravity measuring devices,

d) solving said equations to derive values representative of all of said rotation angles, using measured values of said outputs from said first and second sets of gravity component measuring devices, and

e) computing from said values of said rotation angles the differential azimuth angle between said first set and said second sets of gravity measuring devices.

A first alternative embodiment is shown that is computationally simpler but has certain geometrical errors for larger bend angles. These geometrical errors have a well defined characteristic and corrections may be computed for some applications.

A second alternative embodiment is also shown that is even simpler computationally than either of the first two embodiments. This alternative has geometric errors comparable to those of the first alternative embodiment but the errors do not have a well defined characteristic and are not suitable to computed corrections.

Another object is to provide in either of such methods for inclusion of the first set of gravity measurement devices with a measure-while-drilling (MWD) apparatus in a drill string, and said second set of said gravity component measurement devices is located near the bit of said drill string, and the said computed differential azimuth angle between said first and second gravity component measuring devices is added to the azimuth data from said MWD apparatus to determine the total azimuth angle at the location of said second set near the bit.

A further object is to derive bending in the drill string between the MWD apparatus and the drill bit from the described total azimuth angle.

DRAWING DESCRIPTION

FIG. 1 shows apparatus for accurately deriving differential azimuth, according to the invention;

FIG. 1 a shows FIG. 1 type apparatus in an MWD environment;

FIG. 2 illustrates two sets of gravity component measurement devices, in a borehole;

FIG. 3 shows prior apparatus.

DETAILED DESCRIPTION

U.S. Patent Application Ser. No. 06/200,096, “Borehole Survey Apparatus and Method”, filed Oct. 23, 1980 discloses a concept of using two sets of accelerometers spaced apart some distance along the borehole axis as by intervening structure. The structures between the two sets of accelerometers is flexible to follow changes in direction of the borehole about axes normal to the borehole but is torsionally rigid about the borehole axis, so that there is no change in the relative orientation of the two sets about the borehole axis.

A description of the azimuth angle computation is found in that application between page 9, line 29 and page 21, line 25. This long section describes the mechanization used to determine inclination and azimuth differences from the sensed gravity components from the two spaced-apart sets of accelerometers. The method is described in much more detail. The first result is an update to a matrix relating the attitude of the borehole to an earth-fixed coordinate system. See Page 19, Line 14, for an equation for azimuth, based on elements of the matrix, which is believed to be correct. There is an apparent error in the mechanization, which may be a typographical error. At Page 12, Line 25, the value of 7 is defined as:

-   -   γ=α−β         At Page 15, Line 21, the value of γ is defined as:     -   γ=α−δ

The description of the mechanization is quite extensive and entirely different from the very limited disclosure of U.S. Pat. No. 6,480,119.

U.S. Pat. No. 6,480,119, “Surveying a Subterranean Borehole Using Accelerometers”, describes and claims a survey method using two spaced-apart acceleration sensor means to determine change in azimuth between the two sets as the combined sets are traversed along the borehole. FIG. 3 herein is a figure from that patent.

The following description is found in that patent (there appear to be some typographical problems in the U.S. patent, as issued, between the claims and the text).

-   -   “With respect to the present invention, the reference point can         be taken at the surface or at any point on the measurement line.         By way of example we will consider the reference to be from the         surface. From the reference point the known azimuth of the         borehole is considered to be the azimuth of the long axis (the         z-axis) of one of the sets of accelerometers which are         positioned at the same location as the reference, i.e., the top         set. By so considering then the x, y plane and the z axis pole         of the top set of accelerometers now have a fixed spatial         orientation. If the lower set of accelerometers also had the         same orientation then it would give the same sensor outputs, and         no change in either inclination or direction would be seen         between the two sets.”

Now, with reference to FIG. 3, consider what happens if the lower accelerometer set 20 changes direction by an amount shown as ‘delta-azimuth’ in FIG. 3. If the lower accelerometer set 20 had complete freedom with respect to the upper set 30 then it could still produce the same output values and the azimuth change would not be detectable. However, since the two respective accelerometer sets are joined by a structure 10 which allows only bending and not rotation along the long axis, such as a tube, then the lower accelerometer set's x, y plane cannot change dip direction without a change in some of the x, y and z accelerometer values of the set.

The method of determining the change in azimuth of the lower accelerometer set 20, and hence of the borehole, is given by the following relationships.

First, the borehole inclination (Inc1) at the position of the upper accelerometer set 30 must be found. This is given by: Inc1=arc tan(({square root}(Gx1² +Gy1²))/Gz1)  Equation 1

Next, the borehole inclination (Inc2) at the position of the lower accelerometer set 20 must found. Inc2 is given by: Inc2=arc tan(({square root}(Gx2² +Gy2²))/Gz1)  Equation 2

An intermediate value beta is then found equal to: $\begin{matrix} {{Beta} = \frac{\arctan\left\{ {\left( {{{Gx2}*{Gy1}} - {{Gy2}*{Gx1}}} \right)*\left. \sqrt{}\left( {{Gx1}*{Gy1}*{Gz1}} \right) \right.} \right\}}{\left\{ {{{Gz2}\left( {{Gx1}^{2} + {Gy1}^{2}} \right)} + {{Gz1}\left( {{{Gx2}*{Gx1}} + {{Gy2}*{Gy1}}} \right)}} \right\}}} & {{Equation}\quad 33} \end{matrix}$

And having found the above values, the borehole azimuth at the lower accelerometer set 20 can be found by Borehole Azimuth at position of 20=Reference Azimuth+(Beta/(1−Sin(((Inc(1)+Inc(2))/2))))  Equation 4

Having now established the Azimuth and Inclination at the position of the lower accelerometer set (hereafter referred to as position 2), the spatial co-ordinates of the borehole between the position of the upper accelerometer set (hereafter referred to as position 1) and position 2 can be calculated using standard practices e.g. Minimum Curvature calculations. If the top set is moved to a new position between positions 1 and 2 or at position 2 then the new position 1's azimuth becomes the reference Azimuth and the procedure detailed above can be repeated. Note: In Equation 3 the term {square root}(Gx1*Gy1*Gz1) is shown as {square root}(gx1*gy1*gz1) in the patent and application. Since no definition is shown for the lower case g's, then it is assumed that the upper case is intended. The form with the capital letter notation also is shown in a corresponding PCT application. Also, some of the {square root} signs are missing in the patent specification but are included in the claims. No definition is shown for the “intermediate value beta.” Equation 3 is believed to be in error. Those skilled in the arts will recognize at least three problem areas:

1. For an arc tan function, the operand should be dimensionless. In Equation 3 whatever units are used, the result is not dimensionless. This objection could be removed if all of the accelerometer measurements are assumed to be normalized to a local one g field and thus equivalent to direction cosines relating the accelerometer input axes to the local vertical gravity vector.

2. The factor in Equation 3 of gx1*gy1*gz1 is shown in the claim as having a square root. This is believed to be the desired form since the application presented it that way and PCT publication W0 00/11316 which is a corresponding PCT application shows that form. If so, then as the tool rotates about the borehole z axis, either gx1 or gy1 could become negative making the square root a negative number.

3. The same gx1*gy1*gz1 factor can become zero whenever the x or y axis is horizontal in the number 1 position. This would lead to the value of Beta being zero, independent of the value of any measurements at the number two position.

Further, in Equation 4 the term (1−Sin(((Inc(1)+Inc(2))/2)))) is an unusual form to find in an angular-relation problem. It is not unknown. The coversed sine (coversine A=(1−sin A)) is related to the versed sine (ver sine A=(1−cos A)). This later form was used in celestial navigation in the form of the half versed sine (haversin A=½(1−Cos A) in what was called the cosine-haversine formula for computing the elevation of a star above the horizon. It is alleged that this formula was better in the use of logarithm tables and such were published for navigation use.

The mechanizations that follow overcome the problems of the cited prior art. If one considers the Gx1,Gy1 and Gx2,Gy2 pairs to be orthogonal components of the cross-borehole component of the earth's gravity vector, the term Gx2*Gy1−Gy2*Gx1 can be recognized as the cross-product of the two cross-borehole vectors and proportional to the sine of the angle between them. The term Gx2*Gx1+Gy2*Gy1 can be recognized as the dot-product of the two cross-borehole vectors and proportional to the cosine of the angle between them. The arc tan of the ratio of these two quantities would provide the difference in highside angle between the two sets of measurements. The term Gx1²+Gy1² can be recognized as the square of the magnitude of the cross-borehole vector at position 1. No geometrical significance can be seen in the factor {square root}(Gx1*Gy1*Gz1).

Current Invention

Preferred Embodiment

Since the method of this invention determines only changes in the azimuth orientation of a borehole between two spaced-apart locations along the borehole, the initial azimuth angle may be considered as zero. This simplifies the equations below without any loss in the general usefulness of the method.

Given an X,Y,Z earth fixed coordinate set having X horizontal at any arbitrary, assumed zero, azimuth angle with respect to north, Y horizontal and normal to X and Z down and given an x₁,y₁,z₁ borehole fixed coordinate set for a first set of gravity component measurement devices related to X,Y,Z by rotations:

Alpha, α, about Y which is a tilt or inclination angle

Beta, β, about the new Z′ axis which is a gravity highside angle

The x₁,y₁,z₁ coordinate set is related to the X,Y,Z coordinate set by the direction cosine matrix C₁ whose elements are:

-   -   C₁(x₁,X)=Cos(α)Cos(β)     -   C₁(x₁,Y)=Sin(β)     -   C₁(x₁,Z)=−Sin(α)Cos(β)     -   C₁(y₁,X)=−Cos(α)Sin(β)     -   C₁(y₁,Y)=Cos(β)     -   C₁(y₁,Z)=Sin(α)Sin(β)     -   C₁(z₁,X)=Sin(α)     -   C₁(z₁,Y)=0     -   C₁(z₁,Z)=Cos(α)         And given another x₂,y₂,z₂ borehole fixed coordinate set for a         second set of gravity component measurement devices related to         the prior x₁,y₁,z₁ borehole fixed set by rotations:     -   Tau, τ, about the z1 axis which is a bend direction angle of the         drill string     -   Kappa, κ, about the new x₁′ axis which is a bend magnitude angle         of the drill string orthogonal to tau     -   Tau, −τ, about the new z₁″ axis, such rotation being the         negative of the rotation tau above         The x₂,y₂,z₂ coordinate set is related to the x₁,y₁,z₁         coordinate set by the direction cosine matrix C₂ whose elements         are:     -   C₂(x₂,x₁)=Cos(κ)+(1−Cos(κ))Cos²(τ)     -   C₂(x₂,y₁)=(1−Cos(κ))Sin(τ)Cos(τ)     -   C₂(x₂,z1)=−Sin(τ)Sin(κ)     -   C₂(y₂,x₁)=(1−Cos(κ))Sin(τ)Cos(τ)     -   C₂(y₂,y₁)=Cos(κ)+(1−Cos(κ))Sin²(τ)     -   C₂(y₂,z₁)=Cos(τ)Sin(κ)     -   C₂(z₂,x₁)=Sin(τ)Sin(κ)     -   C₂(z₂,y₁)=−Cos(τ)Sin(κ)     -   C₂(z₂,z₁)=Cos(κ)         The sequence of rotation, tau, kappa, −tau, was chosen to         represent a rotation sequence of a drill string such that there         was no rotation about the z, axis between the x₁,y₁.z₁ set of         axes and the x₂,y₂,z₂ set of axes and that the net rotation         angle of the drill string was represented by kappa. For a         direction cosine matrix such as C₂ which represents a rotation         of a set of coordinates, it may be shown that the z-axis         component of the net rotation is given by C₂(y₂,x₁)−C₂(x₂,y₁)         which can be seen to be zero for the values shown above. Further         for such a direction cosine matrix, C₂, it may be shown that         cosine of the total rotation angle, in this case the intended         kappa, is given by:     -   Rotation Angle=arc cos(Tr(C₂)−1)/2)         where Tr(C₂), the trace of the matrix is given by         C₂(x₂,x₁)+C₂(y₂,y₁)+C₂(z₂.z₁). Using this relation, the Rotation         Angle becomes the arc cos(Cos(κ)) which is just kappa, κ, as         desired.

Given these relationships, consider two orthogonal sets of accelerometers, spaced apart along the z axis. The bend direction angle tau and the bend magnitude angle kappa represent the rotations of the lower set with respect to the upper set. There is no rotation about the borehole z axis between the two measurement sets since they are connected by a structure which is torsionally rigid about the borehole axis. Consider the upper set as Set 1 and the lower set as Set 2.

For Set 1, the inclination angle is alpha and the highside angle is beta. Consider the magnitude of gravity to be one unit. The three gravity components of Set 1 then are:

-   -   G_(x1)=−Sin(α)Cos(β)     -   G_(y1)=Sin(α)Sin(β)     -   G_(z1)=Cos(α)         From these, one may compute α and β as:     -   α=arc tan(({square root}G_(x1) ²+G_(y1) ²)/G_(z1)) and β=arc tan         (G_(y1)/−G_(x1))         For Set 2, the gravity components are measured. Their values can         be expressed in terms of a second set of rotations, equivalent         to alpha and beta for Set 1 and a tilt or inclination value, α₂,         and a gravity highside angle, β₂, can be computed from the         measurements. What is needed to solve for a change in azimuth         angle between Set 1 and Set 2 is the direction cosine matrix C₃         which relates the x₂,y₂,z₂ coordinate set to the original         starting coordinate set X,Y,Z.         The direction cosine matrix C₃ is by definition equal to the         matrix product of the matrix C₁ times C₂ so since the matrix C₁         is known from the Set 1 accelerometer measurements the problem         becomes that of finding the matrix C₂ from the Set 1 and Set 2         gravity component measurements.         The gravity components of G₂ for Set 2 are related to the         gravity components G₁ for Set 1 by the matrix relationship:     -   [G₂]=C₂[G₁]         where [G₂] has three components Gx₂,Gy₂ and Gz₂, C₂ is the         direction cosine matrix defined above relating the x₂,y₂,z₂         coordinate set and [G₁] has three components G_(x1),G_(y1), and         G_(z1).         The problem is, given the measurements [G₁], [G₂], and the form         of C₂ in terms of the defined angles tau and kappa, solve for         numerical values of tau and kappa. From these, the required         values in the matrix C₂ can be found and then the required         values from the matrix C₃ can be found.         For convenience, define the vector [d₁] as equal to [G₁] and         [d₂] as equal to [G₂].         Further define the vector [n₁] as equal to its three components         [C₁(x₁,X),C₁[y₁,X],C₁[z₁,X] and further define the vector [e₁]         as equal to its components [C₁(x₁,Y),C₁(y₁,Y),C₁(z₁,Y)]. Thus         the vectors [n₁], [e₁] and [d₁] are just the X, Y and Z axes of         the original starting point coordinates expressed in the         x₁,y₁,z₁ coordinate set.         The vector dot product [d₁]·[d₂], or equivalently [G]·[G₂], the         gravity measurement vectors, is equivalent to [d₁]·[C₂[d₁]].         Since the form of C₂ is known in terms of the angles tau and         kappa the following equations related to [d₁]·[d₂] can be         written using alpha, beta, tau and kappa as defined above. The         values for alpha and beta are known from the Set 1 measurements         and tau and kappa are the unknowns for which a solution is         desired.     -   [d₁]·[d₂]=[d₁]·[C₂[d₁]]=Cos(κ)+(1−Cos(κ))Sin²(α)Cos²(β+τ)         from which     -   Cos(κ)=([d₁]·[d₂]−Sin²(α)Cos²(β+τ))/(1−Sin²(α)Cos²(β+τ))         Similarly,     -   [n₁]·[d₂]=[n₁]·[C₂[d₁]] and [e₁]·[d₂]=[e₁]·[C₂[d₁]]         Thus,     -   [n₁]·[d₂]=(1−Cos(κ))Sin(α)Cos(β+τ)){−Cos(α)Cos(β+τ)}−Sin(κ)Sin(β+τ)     -   [e₁]·[d₂]=−(1−Cos(κ))Sin(α)Cos(β+τ)){         Sin(β+τ)}−Sin(κ)Cos(α)Cos(β+τ)         To simplify this, define     -   p=−Cos(α)Cos(β+τ) and q=Sin(β+τ)         then     -   [n₁]·[d₂]=(1−Cos(κ))Sin(α)Cos(β+τ))p−(Sin(κ))q         and     -   [e₁]·[d₂]=−(1−Cos(κ))Sin(α)Cos(β+τ))q−(Sin(κ))p         since     -   p²+q²=1−(1−Sin²(α)Cos²(β+τ))         then     -   ([n₁]·[d₂])q+([e₁]·[d₂])p=−(Sin(κ))(p²+q²)         and     -   Sin(κ)=−(([n₁]·[d₂])q+([e₁]·[d₂])p)/((1−Sin²(α)Cos²(β+τ))         From Sin(κ) and Cos(κ) the value of kappa, κ, can be computed.         A similar treatment could be used to find a term including tau         by eliminating all terms in kappa from the above equations.         However, a simpler solution is found by substituting for both         Sin(κ) and Cos(κ) into the equation for [n₁]·[d₂] shown above.         These substitutions result in an equation such that     -   ([n₁]·[d₂])p−([e₁]·[d₂])q=(1−[d₁]·[d₂])Sin(α)Cos(β+τ))         Using the definitions w_(n)=([n₁]·[d₂]), w_(e)=([e₁]·[d₂]) and         w_(d)=[d₁]·[d₂] and expanding out the equations for p and q         shown above, the result is     -   w_(e) Sin(β+τ)=[Cos(α)w_(n)+(1−w_(d))Sin(α)] Cos(β+τ)         For this equation to be true, the coefficient of Sin(β+τ) must         be proportional to Cos(β+τ) and vice versa:     -   Sin(β+τ)=K[Cos(α)w_(n)+(1−w_(d))Sin(α)] Cos(β+τ)=−K w_(e)         where K is some proportionality constant. From these relations,         since α and β are known from the Set 1 measurements one can         solve for τ and then given this result, κ can be found using the         equations above for Sin(κ) and Cos(κ). Since K can be positive         or negative, there is an ambiguity as to the proper value for τ.         From physical reasoning, the value of τ that results in the         smallest κ should be chosen. Alternatively, returning to the         equation shown above:     -   w_(e) Sin(β+τ)=[Cos(α)w_(n)+(1−w_(d))Sin(α)] Cos(β+τ)         This can be rewritten as:     -   A Cos(β+τ)+B Sin(β+τ)=0         which after dividing through by A²+B² is equivalent to     -   Sin(θ+β+τ)=0         where θ=arc tan(A/B) and thus τ=−(θ+β) and again their may be an         ambiguity since Sin(180)=Sin(0). The alternative expression         should be τ=180−(θ+β). As before, the value of tau that leads to         the smallest kappa should be used. Having found values for kappa         and tau using the measured data from the Set 1 and Set 2         accelerometers, the values for the nine elements of the         direction cosine matrix C₂ may be computed. From this, and from         the values for the direction cosine matrix C, determined from         values for alpha and beta derived from the Set 1 accelerometer         measurements, the matrix product:     -   C₃=C₂C₁         may be found. From the direction cosine matrix C₃, a result for         the azimuth orientation of the Set 2 measurements may be found.         Note that by definition, the azimuth direction of a borehole is         the horizontal projection of the along-borehole z axis onto the         horizontal plane defined by the X,Y, vectors cited above. Given         the values for the direction cosine matrix C₃, the azimuth angle         change, ΔΨ, from Set 1 to Set 2 is then found as:     -   ΔΨ=arc tan(C₃(z₂,Y)/C₃(z₂,X))         If ΔΨ is added to some initial azimuth angle for Set 1, then the         result is the total azimuth angle for Set 2.

First Alternative Embodiment

Given an X,Y,Z earth fixed coordinate set having X at any arbitrary angle with respect to north, Y normal to X and Z down and given an x,y,z borehole fixed coordinate set for a first set of gravity component measurement devices related to X,Y,Z by rotations:

-   -   Alpha, α, about Y which is a tilt angle     -   Beta, β, about the new Z′ axis which is a highside angle         And given another x,y,z borehole fixed coordinate set for a         second set of gravity component measurement devices related to         the prior x,y,z borehole fixed set by rotations:     -   Gamma, γ about the new Y″ axis which is a bend angle of the         drill string     -   Delta, δ, about the new X′″ axis which is a bend angle of the         drill string orthogonal to gamma         This last x,y,z coordinate set is related to the X,Y,Z         coordinate set by a nine-element direction cosine matrix. The         elements of that matrix are:

-   C(x,X)=−Sin(α)Sin(γ)+Cos(α)Cos(β)Cos(γ)

-   C(x,Y)=Sin(β)Cos(γ)

-   C(x,Z)=−(Cos(α)Sin(γ)+Sin(α)Cos(β)Cos(γ))

-   C(y,X)=Sin(δ)(Sin(α)Cos(γ)+Cos(α)Cos(β)Sin(γ))−Cos(α)Sin(β)Cos(δ)

-   C(y,Y)=Cos(β)Cos(δ)+Sin(β)Sin(γ)Sin(δ)

-   C(y,Z)=Sin(α)Sin(β)Cos(δ)+Sin(δ)(Cos(α)Cos(γ)−Sin(α)Cos(β)Sin(γ))

-   C(z,X)=Cos(α)Sin(β)Sin(δ)+Cos(δ)(Sin(α)Cos(γ)+Cos(α)Cos(β)Sin(γ))

-   C(z,Y)=Sin(β)Sin(γ)Cos(δ)−Cos(β)Sin(δ)

-   C(z,Z)=Cos(δ)(Cos(α)Cos(γ)−Sin(α)Cos(β)Sin(γ))−Sin(α)Sin(β)Sin(δ)     Given these relationships, consider two orthogonal sets of     accelerometers, spaced apart along the z axis. The bend angles gamma     and delta represent the rotations of the lower set with respect to     the upper set. There is no rotation about the borehole z axis     between the two measurement sets since they are connected by a     structure which is torsionally rigid about the borehole axis.     Consider the upper set as Set 1 and the lower set as Set 2.     For Set 1, gamma and delta are both zero. The inclination angle is     alpha and the highside angle is beta. Consider the magnitude of     gravity to be one unit. The three gravity components of Set 1 then     are:     -   G_(x1)=−Sin(α)Cos(β)     -   G_(y1)=Sin(α)Sin(β)     -   G_(z1)=Cos(α)         From these, one may compute α and β as:

-   α=arc tan(({square root}G_(x1) ²+G_(y1) ²)/G_(z1)) and β=arc     tan(G_(y1)/−G_(x1))     For Set 2, substituting the above gravity components for the like     indicated relations:     -   G_(x2)=−(G_(z1) Sin(γ)−G_(x1) Cos(γ))     -   G_(y2)=G_(y1) Cos(δ)+(G_(x1) Cos(γ)−G_(x1) Sin(γ) )Sin(δ)     -   G_(z2)=(G_(z1) Cos(γ)+G_(x1) Sin(γ))Cos(δ)−G_(y1) Sin(δ)         The immediate problem is, given the two sets of gravity         measurements, what are the values of gamma and delta from the         last three equations above. Since there are really only two         variables, gamma and delta, and there are three equations in         this second set a solution is available.         The first equation is recognized as equivalent in form to         Equation 401.2 in “Tables of Integrals and Other Mathematical         Data,” Herbert Bristol Dwight, The Macmillan Company, New York,         N.Y., 1955. That equation is:     -   p Cos A+q Sin A=r Sin(A+θ)         where     -   r={square root}(p²+q²)Sin θ=p/r Cos θ=q/r         Thus the equation:     -   G_(x2)=−(G_(z1) Sin(γ)−G_(x1) Cos(γ))         may be equated to the form above wherein p=G_(x1), q=−G_(z1),         from which     -   r={square root}(G_(x1) ²+G_(z1) ²) and θ=arc tan(G_(x1)/−G_(z1))         This leads to:     -   γ=arc sin(G_(x2)/{square root}(G_(x1) ²+G_(z1) ²)−arc         tan(G_(x1)/−G_(z1))         Since some computer languages do not include an arc sin routine,         this expression can be shown as a arc tan function:     -   γ=arc tan(G_(x2)/{square root}(G_(x1) ²+G_(z1) ²−G_(x2) ²)−arc         tan(G_(x1)/−G_(z1))         Some care is required in choosing the sign for the square-root         function shown. The simplest approach is to compute a value of         delta with both a positive and with a negative value for the         square root. The smallest value will be the one desired.         Similarly in the equation:     -   G_(y2)=G_(y1) Cos(δ)+(G_(z1) Cos(γ)+G_(x1) Sin(γ))Sin(δ)         the form (G_(z1) Cos(γ)+G_(x1) Sin(γ)) may be assigned a value         G_(c) which may be computed since all of the elements G_(z1),         G_(x1) and γ are known. Thus this equation can be put in the         form above wherein p=G_(y1), q=G_(c), r={square root}(G_(y1)         ²+G_(c) ²) and θ=arc tan(G_(y1)/G_(c)).         Using these then the equation may be solved as:     -   δ=arc sin(G_(y2)/{square root}(G_(y1) ²+G_(c) ²)−arc         tan(G_(y1)/G_(c))         Since some computer languages do not include an arc sin routine,         this expression can be shown as an arc tan function:     -   δ=arc tan(G_(y2)/{square root}(Gy_(y1) ²+G_(c) ²−G_(y2) ²)−arc         tan(G_(y1)/G_(c))         Some care is required in choosing the sign for the square-root         function shown. The simplest approach is to compute a value of         delta with both a positive and with a negative value for the         square root. The smallest value will be the one desired.         This completes the solution. However for a check on the value of         δ the third equation above:     -   G_(z2)=(G_(z1) Cos(γ)+G_(x1) Sin(γ))Cos(δ)−G_(y1) Sin(δ)         may be solved as above. Note that G_(c) is the same as above.         Now p=G_(c), q=−G_(y1), r={square root}(G_(y1) ²+G_(c) ²) and         θ=arc tan(G_(c)/−G_(y1)).         Using these then:         δ=arc sin(G_(z2)/{square root}(G_(y1) ²+G_(c) ²)−arc         tan(G_(c)/−G_(y1))         Since some computer languages do not include an arc sin routine,         this expression can be shown as a arc tan function:     -   δ=arc tan(G_(z2)/{square root}(G_(y1) ²+G_(c) ²−G_(z2) ²)−arc         tan(G_(y1)/−G_(c))         Again, some care is required in choosing the sign for the         square-root function. Given now α, β, γ and δ the direction         cosines C(z,X) and C(z,Y) may be computed as:

-   C(z,X)=Cos(α)Sin(β)Sin(δ)+Cos(δ)(Sin(α)Cos(γ)+Cos(α)Cos(β)Sin(γ))

-   C(z,Y)=Sin(β)Sin(γ)Cos(δ)−Cos(β)Sin(δ)     The azimuth angle change, ΔΨ, from Set 1 to Set 2 is then found as:     -   ΔΨ=arc tan(C(z,Y)/C(z,X))         If ΔΨ is added to some initial azimuth angle for Set 1, then the         result is the total azimuth angle for Set 2.

Second Alternative Embodiment

An approximate solution to the problem described above is found from U.S. Pat. No. 4,071,959, Gyro Stabilized Single-Axis Platform, at Column 3, Line 5, to Column 5, Line 51. This patent provides the basis for the Sperry Sun B.O.S.S. survey tool. That tool has a single set of accelerometers mounted on a gyro stabilized single axis platform having its stabilization axis along the borehole. As the tool travels along the borehole, it follows the bends of the borehole but is, in effect, torsionally rigid about the borehole axis. Thus subsequent sets of accelerometer measurements are equivalent to the mode described above. This solution is based on the small-angle assumption for which the sine of the angle is the angle and the cosine of the angle is one.

Equation E of the patent shows the solution as:

-   -   ΔΨ=−(1/Cos Θ))ΔΦ₂         where ΔΨ is the change in azimuth between the two sets of data     -   Θ is the inclination angle     -   ΔΦ₂ is the change in high side angle between the two sets of         data         Placed in the context of the notations shown at the beginning of         this paper this is equivalent to:     -   ΔΨ=−(1/Cos(α))arc         tan((G_(y1)G_(x2)−G_(y2)G_(x1))/(G_(x1)G_(x2)+G_(y1)G_(y2)))         Error Consideration

Simple physical reasoning shows that none of the embodiments cited above produce a useful result if the borehole is in an initially horizontal plane and the angular bends are in that plane. This is simply the consequence of the fact that the output of a set of accelerometers is invariant for a rotation about the gravity vector, the nominally vertical direction. Beyond this limitation, there may/are also be geometrical errors due to the mechanization and as a practical matter accelerometer errors that will cause errors in the gravity component measurements. Such accelerometer errors are not discussed in detail but it is apparent that their effect will be largest for a nominally vertical starting orientation of the borehole.

Detailed error analyses of the geometrical errors in the above methods show that for the preferred embodiment, geometrical errors are zero except for the exception for horizontal boreholes cited above. These same analyses show that the geometrical errors for the first alternative embodiment increase as the magnitude of the bend angle between the two sets of measurements increases. Errors are quite small for bend angles of less than about two degrees but increase to a significant level above this value. However, the errors have a definitive pattern and may be shown to be proportional to Sin α Sin β. As such, it is possible to compute a correction that will reduce the error but not eliminate the error.

Further, these error analyses show that the second alternative embodiment has errors comparable to the first alternative embodiment but that the errors do not show a definitive pattern that would suggest correction by computation.

Attempts to analyze the errors in the mechanization shown in U.S. Pat. No. 6,480,119 all failed due to the requirement for the square root of a negative number, namely the product of G_(x1) G_(y1) G_(z1).

It is recognized that in the treatment above, the sequence of the rotations for the first alternative embodiment, gamma and delta, can be reversed without altering the basic concept of the invention. All that is required is to rewrite the equations for solution using the selected order for these rotations.

Applications

The methods of the embodiments shown above can be used for the complete survey of a borehole, given a means to initialize the process to some known azimuth. Such a survey tool eliminates the use of magnetometers or gyroscopes which would reduce the cost, and possibly the size of the survey tool Another method of use would be to provide the upper set of accelerometers as elements of a conventional Measure-While-Drilling (MWD) apparatus. This MWD apparatus would provide an azimuth reference that could be used with the differential azimuath determination from the lower set of accelerometers to obtain the correct total azimuth for the lower set of accelerometers. In a typical drilling procedure, a mud motor or other motor means is often used between the MWD apparatus and a drill bit at the bottom end of the drill string. This often results in a substantial distance along the borehole between the MWD apparatus and the drill bit, leading to some uncertainty with respect to the inclination and azimuth orientation at the bit and in turn uncertainty in the future path of the drill string. Using the methods described above, a second lower set of accelerometers could be installed near the bit and both inclination and azimuth at the bit could be determined.

In all of the embodiments described above, it has been assumed that three orthogonal gravity components are available for the computations described. Since the earth's gravity field is a known constant for any location the required three orthogonal gravity components can be obtained either by a set of three accelerometers or by a set of two orthogonal accelerometers and a third component derived by computation, knowing the two measured components and the know value of the gravity field.

Referring now to FIG. 1, it shows an upper housing 100, a lower housing 101, and a connector 102 connecting 100 and 101. This equipment is suspended by a line or cable 104, in a borehole 103, in the earth. The line travels to a cable winding device 105 at the surface 106, and a power source 105 a is connected to 105 to controllably rotate it, for raising and lowering the equipment 100, 101, and 102, in the borehole as during a survey. An idler 108 turns the cable.

Within 100 is an upper or first set 110 of gravity component measurement devices, as for example three sensitive axis accelerometers, and within 101 is a lower or second set 111 of gravity component measurement devices, as for example accelerometers. Accelerometer sets, generally, are disclosed in Engebretson U.S. Pat. Nos. 5,155,916 and 5,321,893. The accelerometer outputs at both sets are transmitted as via suitable wiring 112 along the cable or line, and to processing circuitry that may include a data recorder 114, and processor 114 for calculating in accordance with the above disclosure. A suitable display 115 and/or keyboard 116 may be connected to 114, for displaying data, or processed data, in accordance with the outputs, including ψ, or Δ Ψ, referred to above, in accordance with the computations referred to, as may be carried out in the processor 114.

Connector 102 is typically such that no net relative rotation occurs as between 110 and 111, abut the borehole axis 103 a, during operation. Borehole engaging guides 120 and 121 on 100 and 101, act to guide 100 and 101 as they move along the borehole, and also to center 100 and 101 in the hole. Guides 120 and 121 may assist in preventing such net relative rotation, as by slipping on the borehole surface. The units 100 and 101 may incorporate suitable power surfaces, if needed. Wiring 112 may provide for up-down control and transmission as between 116 and the equipments 110 and 111.

FIG. 1 a shows the housing 100 associated with or carried by an MWD apparatus 130, on a drill stem 131, in a wellbore 132, and housing 101 associated with or carried by a bit box 144 at the lower end of the stem 131. A section 131 a of stem 131 extends between 100 and 101, and prevents relative rotation as between 100 and 101, as when rotary drilling is interrupted, but the stem being moved axially, to allow measurements. Data from 110 and 111 in 100 and 101 may be transmitted to or processed at the surface or via suitable wiring associated with stem 131, or by other known transmission, for example.

Bending in the drill string between the MWD apparatus and the drill bit is derived during computation of said total azimuth angle. 

1. A method for computing the differential azimuth between two sets of spaced-apart gravity component measurement devices in a borehole, subject to the constraint that there is no net rotation about the borehole axis between the two sets of gravity component measuring devices having outputs, comprising the steps of: a) modeling the relationship of the sensitive axes of a first set of said gravity component measurement devices to an earth fixed coordinate set in terms of an inclination rotation angle, for example alpha, from the vertical and an orthogonal highside rotation angle, for example beta, about said borehole longitudinal axis, b) modeling the relationship of the sensitive axes of a second set of said gravity component measurement devices to said sensitive axes of said first set of gravity measurement devices in terms of a rotation angle, for example tau, a bend direction angle, about said borehole axis, and another rotation angle, for example kappa, a bend magnitude angle, about an axis that is perpendicular to said borehole axis and defined by the bend direction angle tau and a last rotation angle that is the negative of said first rotation angle tau about the said borehole axis, c) forming equations relating said rotation angles alpha, beta, tau, and kappa to the outputs of said gravity measuring devices, d) solving said equations to derive the values representative of all of said rotation angles, using measured values of said outputs from said first and second sets of gravity component measuring devices, and e) computing from said values of said rotation angles the differential azimuth angle between said first set and said second set of gravity measuring devices.
 2. A method for computing the differential azimuth between two sets of spaced-apart gravity component measurement devices in a borehole, subject to the constraint that there is no net rotation about the borehole axis between the two sets of gravity component measuring devices having outputs, comprising the steps of: a) modeling the relationship of the sensitive axes of a first set of said gravity component measurement devices to an earth fixed coordinate set in terms of an inclination rotation angle, for example alpha, from the vertical and an orthogonal highside rotation angle, for example beta, about said borehole longitudinal axis, b) modeling the relationship of the sensitive axes of a second set of said gravity component measurement devices to said sensitive axes of said first set of gravity measurement devices in terms of a rotation angle, for example gamma, about an axis perpendicular to said borehole axis and another rotation angle, for example delta, about another axis that is perpendicular to said borehole axis and also perpendicular to other axis, c) forming equations relating said rotation angles alpha, beta, gamma and delta to the outputs of said gravity measuring devices, d) solving said equations to derive values representative of all of said rotation angles, using measured values of said outputs from said first and second sets of gravity component measuring devices, and e) computing from said values of said rotation angles the differential azimuth angle between said first set and said second sets of gravity measuring devices.
 3. The method of claim 1 wherein said first set of said gravity component measurement devices is included with a measure-while-drilling (MWD) apparatus in a drill string, and said second set of said gravity component measurement devices is located near the bit of said drill string, and the said computed differential azimuth angle between said first and second gravity component measuring devices is added to the azimuth data from said MWD apparatus to determine the total azimuth angle at the location of said second set near the bit.
 4. The method of claim 2 wherein said first set of said gravity component measurement devices is included with a measure-while-drilling (MWD) apparatus in a drill string, and said second set of said gravity component measurement devices is located near the bit of said drill string, and the said computed differential azimuth angle between said first and second gravity component measuring devices is added to the azimuth data from said MWD apparatus to determine the total azimuth angle at the location of said second set near the bit.
 5. A method for computing the differential azimuth between two sets of spaced-apart gravity component measurement devices in a borehole, subject to the constraint that there is no net rotation about the borehole axis between the two sets of gravity component measuring devices having outputs, comprising the steps of: a) obtaining values for an inclination rotation angle α from the vertical, and an orthogonal highside rotation angle β about the longitudinal axis of the borehole, as related to the sensitive axes of a first set of said devices and an earth fixed coordinate set, b) obtaining values for a rotation angle tau, a bend direction angle about the borehole axis, and another rotation angle kappa, a bend magnitude angle about an axis that is perpendicular to the borehole axis, and a last rotation angle that is the negative of tau, as related to the sensitive axes of a second set of said devices. c) deriving values representative of all of said rotation angles, using measured values of said outputs from said first and second sets of gravity component measuring devices, and d) and deriving from said last named values the differential azimuth angle between said first set and said second set of gravity measuring devices.
 6. A method for computing the differential azimuth between two sets of spaced-apart gravity component measurement devices in a borehole, subject to the constraint that there is no rotation about the borehole axis between the two sets of gravity component measuring devices having outputs comprising the steps of: a) obtaining values for an inclination rotation angle α from the vertical, and an orthogonal highside rotation angle β about the longitudinal axis of the borehole, as related to the sensitive axes of a first set of said devices and an earth fixed coordinate set, b) obtaining values for a rotation angle γ about an axis perpendicular to the borehole axis, and another rotation angle Δ about another axis that is perpendicular to the borehole axis and also perpendicular to said other axis, c) deriving values representative of all of said rotation angles, using measured values of said outputs from said first and second sets of gravity component measuring devices, d) and deriving from said last named values of said rotation angles the different azimuth angle between said first set and said second sets of gravity measuring devices.
 7. The method of claim 5 wherein said first set of said gravity component measurement devices is included with a measure-while-drilling (MWD) apparatus in a drill string, and said second set of said gravity component measurement devices is located near the bit of said drill string, and the said computed differential azimuth angle between said first and second gravity component measuring devices is added to the azimuth data from said MWD apparatus to determine the total azimuth angle at the location of said second set near the bit.
 8. The method of claim 6 wherein said first set of said gravity component measurement devices is included with a measure-while-drilling (MWD) apparatus in a drill string, and said second set of said gravity component measurement devices is located near the bit of said drill string, and the said computed differential azimuth angle between said first and second gravity component measuring devices is added to the azimuth data from said MWD apparatus to determine the total azimuth angle at the location of said second set near the bit.
 9. The method of claim 5 wherein said first set of said gravity component measurement devices is part of a measure-while-drilling (MWD) apparatus in a drill string in the borehole and said second set of said gravity component measurement devices is located in the drill string near the bit of said drill string.
 10. The method of claim 9 wherein the said derived differential azimuth angle between said first and second gravity component sensing devices is added to azimuth data derived from said MWD apparatus to determine the total azimuth angle at the location of said second set near the said bit.
 11. The method of claim 10 wherein bending in the drill string between the MWD apparatus and the drill bit is derived from said total azimuth angle.
 12. The method of claim 5 wherein said differential azimuth angle between said sets of devices is further utilized to determine drill string bending between said devices.
 13. The method of claim 1 characterized in that the x₁, y₁, z₁ coordinate set is related to the earth fixed coordinate set X, Y, Z by the direction cosine matrix C₁ whose elements are where: C₁(x₁,X)=Cos(α)Cos(β) C₁(x₁,Y)=Sin(β) C₁(x₁,Z)=−Sin(α)Cos(β) C₁(y₁,X)=−Cos(α)Sin(β) C₁(y₁,Y)=Cos(β) C₁(y₁,Z)=Sin(α)Sin(β) C₁(z₁,X)=Sin(α) C₁(z₁,Y)=0 C₁(z₁,Z)=Cos(α) x₁, y₁, z₁ is a borehole fixed coordinate set for a first set of gravity component measurement devices, alpha, α, about Y, is a tilt or inclination angle, beta, β, about the new Z′ axis, is a gravity highside angle.
 14. The method of claim 2 wherein said differential azimuth ΔΨ is the value arc tan C(z,Y)/C)(z,X) where: C(z,Y)=Sin(β)Sin(γ)Cos(δ)−Cos(β)Sin(δ) and C(z,X)=Cos(α)Sin(β)Sin(δ)+Cos(δ)(Sin(α)Cos(γ)+Cos(α)Cos(β)Sin(γ)) where: α=alpha in claim 2 β=beta in claim 2 γ=gamma in claim 2 δ=delta in claim
 2. 